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ABSTRACT 

A theory  of  system  type  for  linear  mul- 
tivariable servomechanisms  is  developed.  New 
characterizations  cf  system  type  are  formu- 
lated and  two  systematic  methods  for  deter- 
mining system  type  are  derived.  Both  methods 
are  algorithmic  and  easier  to  apply  than  the 
techniques,  given  in  earlier  studies.  The 
first  applies  to  the  case  where  the  open-loop 
system  is  described  by  its  transfer  function 
matrix,  the  second  where  its  time  domain 
state-space  representation  is  given.  New  re- 
sults are  also  presented  for  characterizing 
and  identifying  the  system  type  for  certain 
composite  systems  from  a knowledge  of  the  in- 
dividual subsystem  types. 


I. 


INTRODUCTION 


In  classical  servomechanism  theory  [1, 

2],  "system  type"  is  a well  established  and 
widely  used  concept  for  determining  the 
steady-state  tracking  error  characteristics 
of  a closed-loop  system  from  a knowledge  of 
its  open-loop  transfer  function.  Specifically, 
a single-input/single-output  unity  feedback 
system  (see  Fig.  1)  whose  open-loop  transfer 
function  can  be  put  in  the  form 
KN  ( s ) 

G(s)  = 


s D(s) 


(1) 


is  said  to  be  a "type  t"  servomechanism*.  In 
Fq.  (1),  K is  a nonzero  real  constant;  and 
N(s)  and  D(s)  are  monic  polynomials  in  the 
complex  variable  s which  (i)  are  of  degree  m 
and  n-t,  respectively,  where  si<n  and  0<£<n, 
and  (ii)  have  no  common  factors  or  zero  Foots. 
If  the  closed-loop  system  is  stable,  it  is 
known  [1,2]  that  the  steady-state  error 

e = lim  e(t)  t lim  [y  (t)  - y ( t ) ] 

£-+»  £-»ao  d 

is  zero  for  all  oolynomial  inputs  of  the  form 
f-1  i 

y (t)  = t 

d i=0  i 

where  t^O  and  the  a.  are  arbitrary  real  con- 
stants . 


* It  is  sometimes  also  convenient  [1,2]  to  re- 
fer to  G(s)  as  being  type  l. 


Y (s) 
d 


E (s) 

G (s) 

Y (s) 

f * 

- mJ 

• 

Y - input  (desired  output) 
d 

E - error 


Y - output 
G - transfer 
function 


Fig.  1 Unity  feedback  system. 

The  extension  of  this  concept  to  the 
case  where  the  system  in  Fig.  1 is  a multivari- 
able (m-input/r.-output)  , m>l,  system  is  a non- 
trivial matter  and  has  been  the  subject  of  a 
number  of  investigations  [3-8].  The  difficul- 
ty stems  from  the  interactions  of  the  various 
inputs  with  the  outputs.  This  unfortunately 
precludes  identification  of  system  type  from 
the  open-loop  transfer  function  by  direct  in- 
spection as  in  the  scalar  case.  More  funda- 
mentally, there  is  the  issue  of  what  consti- 
tutes a useful  definition  of  type  for  multi- 
variable  systems. 

Wiberg  [3]  defines  a unity  feedback  mul- 
tivariable servomechanism  to  be  type  t-1  if 
the  m-dimensional  error  vector  e(t)  = y<j(t)- 
y(t)  has  the  property  that 

lim  e (t)  = 0 
when  the  input  is 


y (t)  = 
d 


E-l 


(2) 


where  tiO  and  the  a^  are  arbitrary  real  con- 
stants. He  then  shows  that  a necessary  and 
sufficient  condition  for  the  closed-loop  sys- 
tem to  be  type  i is  that  the  m x m open-loop 
transfer  function  matrix  G(s)  be  expressible 


G(s) 


s 


R (s)  + P(s> 


where  R(s)  and  P(s)  are  such  that 

and  | | lim  s^“* 
s »0 


lim  sR  (s) 


lim 
s - 0 


P(s) 
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whore  (•)“*  and  ||-||  denote  matrix  inverse 
and  norm,  respectively.  Except  when  m is  two 
or  three,  the  above  condition  is  difficult 
and  awkward  to  apply. 

In  : ho  approach  taken  by  Sandell  and 
Athans  ( 4 I {see  also  [5]),  a unity  feedback 
multivariable  servomechan i sn  i-  said  to  be 
type  < if  t is  the  largest  nonnegative  inte- 
ger for  which  the  open-loop  transfer  function 
can  bo  written 

G(s)  - G'  (s) 


where  G' (3)  has  nonsingular  d.c.  gain.  The 
latter  requires  that 


lim 

s»0 


AC'  (s) 


dot  G* (s) 


¥ 0 


where  Af,>  (s)  is  the  characteristic  polynomial 
of  G'(s),  that  is  the  least  common  denomina- 
tor of  all  minors  of  G' (s) , and  det  denotes 
the  determinant.  This  condition  is  equiva- 
lent to  G' (s)  having  no  zeros  [15]  at  the  ori- 
gin. Using  their  definition,  Sandell  and 
Athans  show  that  the  closed-loop  system,  if 
stable,  has  zero  steady-state  error  for  an 
input  given  by  Eq.  (2).  Their  condition  is 
easier  to  check  than  Wiberg's,  but  also  re- 
quires a certain  amount  of  trial  and  error. 

Hosoe  and  I to  [6]  were  the  first  to  ob- 
serve that  the  use  of  a single  number  is  in- 
adequate for  characterizing  system  type  in 
the  multivariable  case.  For  example,  if 


G(s) 


the  closed-loop  system  has  zero  steady-state 
error  when  the  input  to  the  system's  first 
"channel"  is  a step  and  that  to  its  second 
channel  is  a ramp.  As  a result,  an  accurate 
description  of  system  type  in  this  case  is 
[ 8. ! 8 2 1 = [ 1 2]  in  contrast  to  both  the  Wi- 

berg  and  Sandell-Athans  formulations  which 
indicate  that  the  closed-loop  system  is  type 
1.  The  definition  of  system  type  thus  pro- 
posed by  Hosoe  and  Ito  is  in  terms  of  the  vec- 
tor type  [11. ..8  ] where  l^O  is,  for  each 
i=l,2,...,m,  the  system  type  for  the  ith 
input-output  pair  in  the  classical  sense. 

Their  test  for  system  type  consists  of  first 
expressing  G(s)  in  a partial  fraction  expan- 
sion 

P 1 m 

G (s)  = y -4  G + G(s) 

3=1  j 


where  Gj  are  nonzero  constant  matrices,  and 
G(s)  is  rational  and  has  no  poles  at  zero. 
Lettinq  Gl  bo  the  matrix  obtained  by  deleting 
the  it-'1  rdw  of  Gj,  one  forms 


r 

r 

G.  G,  Gp 

Gt  G:  •••  G£ 

N = 

G,  ...  0 

and  N^= 

gJ  g*  . . . 0 

1 

O • 
TZ 

O 

O • « 

_GP  0 •••  0 - 

Then,  * rank  N - rank  Nj . This  test,  which 
is  conceptually  simple,  can  become  involved 
when  p is  large. 

Young  and  Willems  [7]  approach  the  type- 
1 multivariable-  servomechanism  problem  by  aug- 
menting the  usual  time-domain  state-space  des- 
cription (A,U,C)  by  the-  integral  of  the  error, 
viz.  , 


t 

z - / [y  (t ) - y (1 ) Idt 
JO  r 


They  then  show  that  the  augmented  system  is 
controllable  if  and  only  if  (A,B)  is  control- 
lable and 

A B I 


rank 


-C 


= n + m 


With  these  conditions  satisfied,  they  are  able 
to  assign  the  closed-loop  poles  arbitrarily 
via  state  variable  feedback,  and  so  stabiliz- 
ing the  system  guarantee  zero  steady-state 
error  when  yf  is  a vector  of  step  inputs.  The 
extent  to  which  any  channel  of  the  closed-loop 
gives  zero  steadv-st.ate  error  for  higher-order 
inputs,  e.a.,  ramps,  parabolas,  etc.,  is  un- 
known. Porter  and  Bradshaw  [8]  subsequently 
extended  these  results  via  higher  order  aug- 
mentation to  arbitrary  type.  Interestingly, 
their  conditions  for  controllability,  and 
therefore  for  pole  placement,  are  identical  to 
those  above.  Earlier  work  along  these  lines 
had  been  conducted  independently  by  Smith  and 
Davison  [ 9 ] . 

An  alternate  viewpoint  on  the  problem  of 
system  type  is  to  consider  the  system  to  be 
asymptotically  decoupled  for  a certain  class 
of  inputs.  In  this  case,  the  problem  is  simi- 
lar to  those  posed  and  solved  by  Davison  [10] 
and  Huang  [11]  . 

A distinctively  different  approach  to 
the  problem  of  determining  conditions  for  zero 
steady-state  tracking  error  is  given  by  the 
geometric  theory  of  linear  multivariable  con- 
trol [12].  Conditions,  expressed  in  terms  of 
invariant  subspaces,  which  guarantee  the  de- 
sired tracking  properties  have  been  obtained 
by  Bhattacharyya  and  Pearson  [13,14],  Wonham 
[15],  Wonham  and  Pearson  [16],  and  Francis 
and  Wonham  [17],  System  type  is  not  dealt 
with  explicitly  in  these  studies,  but  it 
would  seem  that  an  intimate  connection  should 
exist. 

The  approach  taken  here  is  based  on  a 
generalization  of  the  Hosoe-Ito  definition. 

It  leads  to  new  characterizations  of  system 
type  and  provides  methods  for  determining 
type  from  either  the  open-loop  transfer  func- 
tion matrix  G(s)  or  the  state-space  represen- 
tation (A,B,C,D)  which  are  simpler  to  apply 
than  those  in  previous  works.  These  topics 
are -treated  in  Sections  II  and  III.  In  Sec- 
tion IV,  new  results  are  obtained  for  the 
characterization  and  identification  of  type 
for  composite  systems.  The  conclusion  is  given 
in  Section  V.  In  a companion  paper  [18],  the 
theory  and  methods  for  achieving  a specified 
type  via  cither  pre-  or  post-compensation  are 
presented. 
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TYPE  f1,*, 
MATRICES 


. 1 | TRANSFER  FUNCTION 
m 


The  following  is  a generalization  of  the 
Hosoe-Ito  [6J  definition  of  system  type. 

Definition  1 

A unity  feedback  multivariable  linear 
system  is  called  type  [ 1 1 i 2 ...  f.m)  if  each 
lz,  i-l,2,...,m,  is  the  largest  integer  such 
that,  the  m*m  open-loop  transfer  function  mat- 
rix G(s)  can  be  factored  into 


G(s)  = H (s ) G'  (s) 


where 

H(s)  = diag  [4-  4-  ...  4- 

I s 1 s 2 s m 

and  G’ (s)  has  no  zeros  at  the  origin.  Here, 
G(s)  is  termed  a type  * * * *1^  transfer 

function  matrix.  □ 

This  definition  generalizes  that  in  [6] 
to  the  extent  of  allowing  2^<0.  In  terms  of 
the  steady-state  tracking  properties,  it  is 
of  no  consequence  whether  the  !. | are  less 
than  zero  or  equal  to  zero.  The  magnitude  of 
when  is  less  than  zero,  however,  will 
be  of  importance  in  the  sequel  when  the  pro- 
blem of  system  type  for  series  connections 
of  systems  is  considered. 

It  is  worthwhile  to  note  at  this  point 
that  system  type  for  multivariable  systems 
is  defined  in  terms  of  a left  divisor.  It 
will  be  shown  that  this  leads  to  meaningful 
results  for  steady-state  errors.  The  follow- 
ing simple  example  shows  that  a right  divisor 
will  give  different  results,  in  fact,  results 
which  are  erroneous. 

Example  1.  Let 


where  for  l->0,  the  a^  are  arbitrary  real  con- 
stants and  OSk^cii  with  the  k£  integers;  and 
for  ai=k^=0. 

The  following  definition  and  lemma  will 
be’ useful  in  the  proof  of  the  main  theorem  of 
this  section. 

Definition  2. 

Wifh  G (s) , H (s) , and  G' (s)  as  in  Defini- 
tion 1,  H ( s ) and  H(s)  are  defined  by 

H(s)  = diag  j ■ ■ • ■ n~—  1 

[_  s 1 sl!  s m J 


H (s)  = 


diag  ["  — i — . . . _L_  *1 

Ls*t  Sl+2  Slt>  J 


where  (l,  if  1;  > 0 1 l,  if  t;  < 0 

£ + = and  l~.  = 

i (0  if  L < 0 10  if  ij  1 0 


so  that 


Factoring  G(s)  according  to  the  definition 
yields 


G (s)  = 


i ° 


Thus  G(s)  is  the  open-loop  transfer  function 
matrix  of  a type  [1  0]  system. 

If  now  an  attempt  is  made  to  find  an 
H(s)  which  is  a right  divisor,  none  can  be 
found  except  the  identity.  This  would  indi- 
cate that  G(s)  is  type  [00].  □ 

With  system  type  as  specified  by  Defini- 
tion 1,  it  will  be  shown  below  that  a stable 
servomechanism  with  a type  f 8. , 4 ? ...  4m] 
open-loop  transfer  function  matrix  yields 
zero  steady-state  error  for  input  signals  of 
the  form 


H(s)  = H (s)  H (s)  = H (s)  H (s)  □ 


Note  that  H(s)  has  only, nonnegative  pow- 
ers of  s as  elements,  while  H(s)  involves  only 
nonpositive  powers  of  s.  Now  G(s)  can  be 
factored  as 

a *\j 

G(s)  = H (s)  G'  (s)  = H (s)  G (s) 
where  &(s)  = H(s)  G'(s). 

- 

Lemma  1.  With  G(s)  = H(s)  G(s)  as  above, 
*G  = AH 

Proof . From  Lemma  1 of  [4]  it  follows 
that  Ag  divides  Ah'>q.  Thus  it  suffices  to  show 
that  deg  (AG)  = deg  (A£)  + deg  (AG) . By  the 
definition  of  H(s)  and  ft(s),  each  factor  of  s 
in  Ah  must , correspond  to  factoring  a pole  of 
G(s)  into  H(s).  Otherwise,  AG  would  remain 
constant  while  the  determinant  would  lose  a 
factor  s”1,  adding  a zero  at  s=0,  which  is  not 
allowed.  Thus,  each  factor  s in  Ajj  reduces  the 
order  of  AG  by  one  when  going  from  Ag.  to  form 
AG.  Also,  the  factorization  cannot  introduce 
any  other  new  poles.  Hence, 


dog  (Ag)  = deg  (Afl)  + deg  (Ag) 


I 


A simple  modification  of  the  proof  of 
this  lemma  shows  that  A^.  - A(|  A_,  . 

Now  the  theorem  on  ’steady-state  error 
can  be  stated  and  proved. 

Theorem  1 . 

Let  G(s)  by  a type  U,fj  ...  trans- 

fer function  matrix.  Then  the  stable  unity 
feedback  system  having  G(s)  as  its  open-loop 
transfer  function  matrix  will  track  inputs 
of  tho  form  given  in  Eg.  (3)  with  zero  steady- 
state  error. 


*•-  A„  Ai  del  (11  (s)  ] 


• dettfr1  (s)  + G (s)  ) (9) 


By  construction  of  H(s), 


yields 


Proof . The  Laplace  transform  of  Eq.  (3) 


Y (s)  = _ 
d 


sk,4l  sk,+l 


where  ( ) denotes  the  transpose.  Then  the 
Laplace  transform  of  the  system  error  is 

E(s)  = (I  + G(s)  l"1  Y (s)  (4) 

d 

Now  with  H (s J as  given  in  befinition  2, 
it  is  noted  that 

lim  H” 1 (s)  Y . (s)  = z (5) 


where  zD  is  a finite  m-vector. 

Thus  if  [I  + G(s)]-1  is  factored  as 

d t G (s) ] _ 1 = {H  (s) - (H~ 1 (s)  + G(s)]}-1 


= [H-1  (s)  + G(s)  ] ~ 1 H~  1 (s) 


and  it  can  be  shown  that 

lim[H-1 (s)  + S(s)]*1  = G_ 


where  Gn  has  finite  elements,  then  the  com- 
bination of  Eqs . (4),  (5),  and  (61  can  be  used 
to  show  that 

lim  s E(s)  = lim  s • lim  E(s) 
s-0  s-0  s-0 

= lim  s • lim  [H-1  (s)  +cS(sir‘ 
s-0  s-0 

• lim  (H-'(s)  Yj (s) ] 
s-0 

- lim  s G0  z = 0 

s-0  (7) 

Finally,  Eq.  (7)  can  be  used  in  the  final- 
value  theorem  of  the  Laplace  transform  to 
prove  the  present  theorem. 

First,  Eq.  (6)  must  be  proved.  By  the 
assumed  stability  of  the  unity  feedback  sys- 
tem, the  zeros  of  the  polynomial  A-  det 
(I  + G(s)]  must  be  in  the  open  left-half  com- 
plex plane  (20).  Then 


A„  det  (H(s)J  = 1 (10) 

Combining  Eqs.  (101,(9)  and  (8)  yields 

Jim  A^  det  [H-1 (s)  + &(s)J  ? 0 (11) 

s-0  G 

With  the  product  in  Eq.  (11)  nonzero,  then 
E<j.  (6)  will  be  verified  if  the  product  of 
Ag  with  each  cofactor  of„[H-‘(s)  + fi(s)]  is  a 
polynomial  in  p.  Since,,  H-  (s)  has  no  poles, 
the  poles  of  (H_I (s)  + &(s)J  are  a subset  of, 
the  poles  of  G(s).  It  then  fol lows,  since  Aq 
contains  a factor  for  each, pole  of  6(s),  that 
Ag  times  any  cofactor  of  [H-1  (s)  + &(s))  is  a 
polynomial.  As  a result, 

lim(H-1  (s)  + fe  (s)  1 ~=' 
s-0  J 

A'i.  cofactor  (H~'(s)  + G(s)].. 

= lim  

s-0  Aq  det  [H-1  (s)  + &(s) ] 

= a finite  number  for  all  i,  j 

where  liiim  and  lijfm.  Thus  Eq.  (6)  is  veri- 
fied and  Eq.  (7)  holds.  Using  the  final-value 
theorem,  the  steady-state  error  is  zero  if  the 
limit  exists.  However,  by  the  assumed  stabil- 
ity of  the  unity  feedback  system,  the  limit  does 
exist  and  the  theorem  is  proved.  0 

As  an  immediate  consequence  of  Lemma  1, 
one  has  the  following  proposition  concerning 
bounds  on  the  type  of  a transfer  function  type. 

Proposition  1.  Let  G(s)  be  a type 
[l  l . . . , icT]  transfer  function  matrix.  Also, 


m 


and  L~  *=  £ Then, 

i=l 


(i)  L+  is  less  than  or  equal  to  .the  mul- 
tiplicity of  the  pole  of  G(s)  at  the 
origin 

and 

(ii) L-  is  greater  than  or  equal  to  the 
multiciplicity  of  the  zero  of  G(s) 
at  the  origin. 

Proof.  (i)  If  L+  is  greater  than  the 
multiplicity  of  the  pole  of  G(s)  at  the  origin, 
then  zeros  at  the  origin  have  been  added  to 
G'  (s)  contradicting  the  definition. 

(ii)  If  L“  is  less  than  the 
multiplicity  of  the  zero  of  G(s)  at  the  origin, 
then  a zero  would  remain  at  the  origin  in  G' (s) 
again  contradicting  the  definition.  □ 

III.  IDENTIFICATION  OF  TRANSFER  FUNCTION  TYPE 


lim  A~  det  [I  4 G(s))  i 0 
s-0 


Using  Lemma  1, 


If  the  dimension  m of  the  transfer  func- 
tion matrix  G(s)  is  small,  say  two  or  three, 
the  type  of  the  system  can  be  found  via  trial 
and  error  using  Definition  1.  However,  for 
larger  values  of  m,  the  procedure  would  be 
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extreme  1 y tedious  since  at  each  stage  the 
characteristic  polynomial  must  be  computed. 

In  this  section,  two  simpler  methods  of  com- 
puting the  system  type  ate  given.  The  first 
method  dives  the  system  type  from  a simple 
observation  of  one  of  the  factors  in  a co- 
ptine fact  or  1 z it  ion  of  the  transfer  function 
matrix.  The  other  method  is  an  algorithm 
which  gives  the  system  type  based  on  rank  con- 
ditions using  the  state-space  representation 
(A,B,C,D)  which  is  assumed  to  be  a minimal 
realization.  hose  methods  are  simpler  than 
the  rank  conditions  in  [6],  especially  if 
there  are  large  power  of  s in  the  denominator 
of  any  element  of  G(s). 

Both  methods  considered  give  the  system 
type  only  to  the  point  of  determining  the 


If  a particular  elenent 


is  actually 


negative,  these  methods  will  indicate  that 
the  is. equal  to  zero.  Thus,  only  the  ele- 

ments of  H(s)  are  determined.  This  is,  how- 
ever, the  only  part  which  is  of  concern  when 
the  steady-state  properties  are  considered  as 
was  demonstrated  in  the  proof  of  Theorem  1 . 

The  first  method  is  presented  by  way  of 
the  following  theorem. 

Theorem  2. 

Let  G(s)  be  the  open-loop  transfer  func- 
tion matrix  of  a unity  feedback  system  and 
let  D-,(s)N(s)  be  a coprime  factorization  of 
G(s).  Then  the  unity  feedback  system  is  type 
[2*1+  ...  2+1  where  the  2t,  i = 1,2,  ...,  m, 
ar£  the  greatest  power  of  s such  that  each 
element  of  column  i of  D(s)  is  divisible  by 
s. 

Proof.  From  Definition  1,  each  2t  is 
naxj.mal  such  that  no  zeros  are  introduce^ 
in  G(s)  by  the  factorization  G(s)  = H(s)G(s). 
Further,  G(s)  can  be  expressed  in  the  form 

G (s)  = H” 1 (s)  G (s) 

= ft"1 (s)  D"1 (s)  N (s) 

where  D_l(s)N(s)  is  any  coprime  factorization 
of  G(s),  and  fi-l(s),  N(s)  and  D(s)  are  all 
iTi'm  polynomial  matrices.  The  polynomial  sys- 
tem matrix  for  G(s)  is 


P(s)  = 


D (s)  N(s) 
H-1 (s)  0 


The  zeros  of  G(s)  are  the  zeros  of  det 
P(s)  after  the  decoupling  zeros  [19]  have 
been  removed.  Since  D(s)  and  N(s)  are  co- 
prime, the  only  decoupling  zeros  are  output 
decoupling  zeros.  Let 

D (s)  N (s) 

P (s)  = 

1 H7‘(s)  0 

denote  P(s)  with  its  output  decoupling  zeros 
removed.  . 

The  zeros  of  G(s)  are  the  solutions  of 
0 = det  P (s) 

i 

« dot  fr 1 (s)  • det  N (s) 

So  in  order  that  no  zeros  be  introduced 


into  G(s),  all  of  the  factors  of  s in  H'‘(s) 
must  be  output  decoupling  zeros  and  hence  re- 
movable. By  repeated  application  of  the  lem- 
ma in  the  Appendix  to  (b+(s)  j il~ 1 (s)T)  with 
'i-0,  the  maximal  powers  of  s which  can  be  re- 
moved from  P(s)  are  the  largest  powers  of  s, 
say  ij,  which  divide  each  element  of  column  j 
of  D ( s j . So  the  largest  power  of  s in  entry 
i of  H-1 (s)  such  that  no  zeros  are  introduced 
in  G(s)  and  the  largest  power  of  s which  ai- 
vides  each  element  of  column  i of  D(u)  are 
the  same,  and  the  theorem  is  proved.  □ 

Algorithms  for  determining  a coprime  fac- 
torization of  a rational  transfer  function 
matrix  can  be  found  in  [22] . The  following 
example  demonstrates  the  method  of  determin- 
ing system  type  which  is  described  in  the 
above  theorem. 

Kxamnle  2.  Let 


G(s)  = 


s2 (s+1) 


A coprime  factorization  yields 


G (s)  = 


(s+1) (s+2 ) 


s2 (s+1) 


s+2  s+2  s 


Then  D(s)  is 


s (s+2) 


0 s ( s+ 1 ) 


(s+1) (s+2) 


s^ (s+1) 


and  it  follows  that  G(s)  is  type  [1  0 2].  □ 

The  second  method  of  determining  transfer 
function  matrix  type  is  for  the  state-space 
representation  (A,B,C,D)  which 'is  assumed  to 
be  a minimal  realization  with  the  four  matrices 
n^m,  n»r,  m*n,  and  m»r,  respectively.  Rather 
than  a direct  characterization  of  type,  the 
following  theorem  presents  an  algorithm  which 
determines  the  value  of  2^+  for  each  i. 

Thcoren  3. 

If  the  state-space  quadruple  (A,B,C,D)  re- 
presents the  open-loop  elements  of  a stable  m- 
input/m-output  unity  feedback  system,  then  the 
following  algorithm  will  determine  the  system 
type . 


i 


( Algorithm.  Let  Q = [CT|  ATCT  | A?TCT! . . . 
,An-lTcT]T(  the  observability  matrix, 

and  let  -..r 

. T * T .Tl  T ,TI  I „ , T .T 
Ql  = • A C1  I A-'  CJ  « An_1  cl 

i t • 3 t iii  i 

Li  , • I J 

where  C ] is  the  matrix  C with  row  i replaced 
by  row  xi  of  CAJ . With  these  matrices,  per- 
form the  following  steps: 

(1)  I.ot  p 5 rank  Q 

(2)  Sot  row  counter  i = 0 

(3)  Increment  row  counter  i by  1 

(4)  Set  j - 0 

(9)  Increment  j by  1 

(6)  Form  0^ 

(7)  Set  k*l  rank  Q* 

(8)  If  k + j = p,  go  to  (5) 

(9)  Set  lt  = j - 1 

(10)  If  i < m,  go  to  (3) 

(11)  End. 


power  of  s.  Then  if  the  rank  Ql  is  decreased 
by  1 when  j is  increased  by  1,  in  output  de- 
coupling zero  was  added.  The  algorithm  thus 
finds  the  largest  j such  that  the  rank  of 
decreases  at  each  step  which,  as  was  seen  u, 
the  proof  of  Theorem  2,  is  the  of  the 
matrix  transfer  function. 

As  was  also  noted  in  the  proof  of  Theor- 
em 2,  the  number  of  output  decoupling  zeros 
introduced  by  element  i of  -H-1  (s)  is  inde- 
pendent of  the  others.  This  justifies  the 
procedure  in  the  algorithm  of  looking  at  one 
row  at  a time.  H 

The  following  example  demonstrates  the 
algorithm. 


Example  3.  Let  the  open-loop  elements 
of  a unity  feedback  system  be  defined  by 
the  quadruple 


r o i o -i 


n o-i 


Proof.  In  the  proof  of  Theorem  2,  it 
was  shown  that  the  tj  of  the  system  type  vec- 
tor are  the  largest  powers  of  s such  that 
multiplication  of  row  i of  G(s)  by  s intro- 
duced exactly  1 1:  output  decoupling  zeros. 
Multiplication  by  larger  powers  would  not  in- 
troduce any  more  decoupling  zeros.  The  ef- 
fect of  multiplication  of  row  i of  G(s)  by  s 
to  some  power  when  G(s)  is  given  by 


A 


0 0 0 


B = 1 0 


Lo  o o 


U i 


0 

0 


0 ■ 

1 . 


0 

0 


Using  the  algorithm,  one  finds  that 


G (s ) = C (si  - A)_IB  + D 


rank  0=3 


is  developed  as  follows. 

The  result  of  multiplying  row  i of  G(s) 
by  si  is  giver,  by  row  i of 

C (si  - A)  ~ 1 B + siD  (12) 

By  using  the  identity  (si  - A) (si  - A)”1 
= I,  the  result 

s (si  - A) ~ 1 = A (si  - A)-1  + I (13) 

is  obtained.  Repeated  application  of  Eq. 

(13)  to  Eq.  (12)  yields  the  result 


rank 

Q.1 

= 

2 

rank 

Qi 

= 

1 

rank 

Q? 

s 

1 

(hence  1 

= 2) 

rank 

Qi 

= 

2 

rank 

Qi 

= 

2 

(hence 

= 1) 

Thus,  [A , B ,C , D]  represents  a type  (2  1]  trans 
fer  function  matrix.  This  can  easily  be  veri- 
fied in  this  case  since 


si  C (si  - A)"1  B + S^D  = CA^ (si  - A) " 1 B 

.k 


. 1-1 
+ T CA  Bs 
k=0 


1 -k-l  j 
J + S D 


G (s) 


The  system  matrix  for  G(s)  with  row  i 
multiplied  by  s3  is  then 


s+1 


□ 


IV.  SYSTEM  TYPE ’FOR  COMPOSITE  SYSTEMS 


Pj  (S) 
1 


Si -A 

ci 


B 

D-i  (S) 
X 


where  C-?  is  given  by  the  algorithm  and  D?  (s) 
is  D with  row  i replaced  by. row  i of  1 


s3  D + s3-k-1  CAkB 

k=0 

From  Theorem  8.1,  Chapter  2 of  Rosen- 
brock  119),  the  rank  defect  of  the  observ- 
ability matrix  Q is  the  number  of  output 
decoupling  zeros.  Thun,  when  j is  increment- 
ed by  1 in  the  algorithm,  this  corresponds  to 
multiplication  of  row  i of  G(s)  by  one  higher 


In  this  section,  some  results  are  obtain- 
ed for  system  type  of  composite  systems.  The 
first  result  is  a bound  on  the  system  typo  for 
the  series  connection  (see  Fig.  2(a)]  in  terms 
of  the  transfer  function  type  of  the  indivi- 
dual transfer  function  matrices.  The  second 
result  gives  conditions  under  which  a trans- 
fer function  type  is  maintained  when  it  is 
connected  in  parallel  with  a constant  gain 
matrix  (see  Fig.  2(b)). 


(a)  Scries  connection 
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1 


(b)  Parallel  connection 
Fig.  2.  Composite  system  connections. 

In  the  first  instance,  one  has  the  fol- 
lowing resalt. 

Theorem  4. 

If  C,  (s)  is  type  [t,Z  ...  *ml  and  G2 (s) 
is  type  {/.  1 1 J . ..tjJ.1  then  the  series  combin- 
ation G,(s)Gj(s)  is  at  least  type  r I , J Z2-> 

. . where  l = min  [t!i;  . (By  at 

least  type  {i,  + t S.?  + Z ...  in+i) , it  is  meant 
that  each  i ? z ^ + f where  tv  is  the  type  for 
the  ith  input-output  pair  of  the  composite 
system.)  The  proof  follows  from  Def.  1.  O 

The  following  corollary  is  an  immediate 
consequence  of  the  theorem. 

Corollary  1.  If  G(s)  is  type  [ 1. 2 fi.  2 ... 
tm)  and  K(s)  is  a dynamic  compensator  which 
is  used  to  stabilize  the  unity  feedback  sys- 
tem with  G(s)*K(s)  as  its  open-loop  transfer 
function  matrix,  then  the  composite  system 
is  at  least  type  [ijij  ...ln)  provided  that 
K(s)  has  no  zeros  at  the  origin. 

Remark.  Corollary  1 was  stated  in  (6J 
without  the  condition  that  K(s)  have  no  zeros 
at  the  origin.  A simple  example  in  the  single 
input/single-output  case  shows  the  necessity 
for  this  requirement. 

The  second  result  of  this  section  per- 


tains to  a type  [!.,£. 


transfer  func- 


tion matrix  which  has  a constant  gain  matrix 
in  parallel  (see  Fig.  2(b)).  This  result  is 
used  in  [23J  in  a study  of  stochastic  multi- 
variable  servomechanisms. 

Theorem  5 ■ 

Let  G ( s ) = D- 1 ( s ) N ( s ) be  a type  ( i , t 2 . . . 
Zml  transfer  function  matrix  and  K an  m*m 
constant  gain  matrix.  Then  G(s)+K  is  at 
least  type  I i , 4 2 ...  Zm)  if  and  only  if  det 
(N(s)+D(s)KJ  has  no  zeros  at  the  origin. 

Proof.  Since  K has  no  poles,  the  poles 
of  the  parallel  combination  G(s)+K  are  simply 
the  poles  of  G(s).  Now  G(s)+K  can  be  factor- 
ed as 


G(s)  + K = H (s ) [G ‘ (s ) + H“*(s)K)  (16) 

and  Eq.  (16)  can  be  written  as 

G(s)  + K = H(s)(G’(s)  + H” 1 .(s)  K ) 

= H(s)D-1 (s) [N  (s ) +D  (s ) H“ 1 (s)K) 

= II  (s)  D~  1 (s)  [N(s)+D(s)K] 

Hence,  G(s)+K  is  at  least  type  [Z,lj  . ..Z  ] 
if  and  only  if  det  (N(s)+D(s)K)  has  no  zeros 


at  r.  = 0.  □ 

Corollary  2.  If  G(u)  is  type  [ 4 , 4 a ... 
Zm)  and  ~Z-  > 0 for  every  i,  then  G(s)tK  is 
at  least  type  (Z,ij  . . . 4m ) . 

Proof.  Using  the  property  of  D(s)  given 
in  Theorem  2,  and  the  assumed  values  for  the 
( every  element  of  the  product  D(s)K  is 
divisible  by  s.  Then  the  constant  term  in 
the  equation  det  (M (s) +D (s) K) =0  must  be  the 
same  as  the  constant  term  in  the  equation 
det  N(s)=0.  Since  N(s)  has  no  zeros  at  the 
origin,  the  constant  term  is  nonzero.  There- 
fore det  [N  (s) +1) (s)  K]  has  no  zeros  at  the  ori- 
gin. □ 

V.  CONCLUSION 

In  this  paper,  a theory  of  system  type 
for  linear  multivariable  servomechanisms, 
which  is  based  on  a generalization  of  the 
Hosoe-Ito  (6)  definition  of  system  type,  has 
been  presented.  The  more  general  definition 
permits  the  development  of  results  for  the 
determination  of  system  type  for  composite 
systems . 

New  characterizations  of  system  type  have 
been  formulated  and  two  systematic  methods 
for  determining  system  type  have  been  derived. 
The  first  method  is  directly  suited  for  appli- 
cation to  the  transfer  function  matrix  (fre- 
quency domain)  description  of  the  system,  the 
second  for  the  state-space  (time-domain)  re- 
presentation. Both  are  easier  to  apply  than 
those  in  previous  studies  and  are  algorithmic 
in  nature.  The  second  method  is  the  only 
known  result  which  treats  multivariable  sys- 
tem type  directly  in  terms  of  the  state-space 
quadruple  (A,B,C,D)  without  need  to  compute 
first  the  corresponding  transfer  function 
matrix. 

The  results  here  on  system  type  for  com- 
posite multivariable  linear  systems  in  terms 
of  the  type  for  the  individual  subsystems  are 
also  new,  and  have  found  applications  else- 
where (18,23). 

APPENDIX 

It  is  known  (19)  that  the  system  matrix 


' T (s)  U(s)  ' 

. -V(s)  W(s)_ 


gives  rise  to  the  transfer  function  matrix 

G(s)  = V(s)  T“‘(s)U(s)  + W(s) 

and  conversely  to  within  strict  system  equi- 
valence where  T,U,V,  and  W are  m^n,  n*r,  m'n, 
and  m*r,  respectively. 

A system  matrix  which  contains  decoupl- 
ing zeros  can  be  reduced  to  one  of  least  or- 
der by  removal  of  the  decoupling  zeros.  The 
following  lemma  concerning  reduction  of  order 
is  needed  in  the  proof  of  Theorem  2. 

Lemma.  Let  U(s)  = diaq  (1  1...1  (s-a)  1 

...  l)  where  (s-a)  is  the  i^h  diagonal  element, 
and  a is  any  complex  number.  Then  rank 
(T(a)  U ( a ) ] < m if  and  only  if  every  element 


I 


w 


of  i of  T(s)  is  divisible  by  s - a. 

Proof.  Sufficiency  is  obvious.  Necess- 
ity comes  f r<  ■:n  the  fact  that  if  any  element 
of  row  i is  not  divisible  by  s - n,  then  at 
s - n the  remaining  m - 1 columns  of  U(u) 
and  one  column  of  T(t)  which  has  a nonzero 
element  in  row  i form  a linearly  independent 
set.  Thus,  rank  [T(s)  U(s)]  = m for  all  com- 
plex r.  if  not  all  elements  of  row  i of  T(s) 
are  divisible  by  s - a . [] 
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